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The properties of fluctuation induced interactions like van dcr Waals and 
Casimir-Lifshitz forces are of interest in a plethora of fields ranging from bio- 
physics to nanotechnology. Here we describe a general approach to compute 
these interactions. It is based on a combination of methods from statistical 

S physics and scattering theory. We showcase how it is exquisitely suited to ana- 
lyze a variety of previously unexplored phenomena. Examples are given to show 
_JL how the interplay of geometry and material properties helps to understand and 

control these forces. 
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1. Introduction 

All material objects, even if charge neutral, support instantaneous current 
fluctuations due to quantum and thermal fluctuations of their charge dis- 
Q\ tribution. The interaction that results from the electromagnetic coupling 

of these currents on different objects is usually called the Casimir force. 
Originally, this force has been derived for two parallel perfect metal plates 1 
and atoms, 2 and generalized later to two infinite dielectric half-spaces with 
planar and parallel surfaces. 3-6 The non-additivity of the Casimir force lim- 
its these results in their applicability to objects at very short separation via 
the so-called proximity force approximation which provides only an uncon- 
trolled approximation of surface curvature to lowest order at vanishingly 
small separations and ignores the global geometrical arrangement of the 
objects. Generically, one encounters in practice geometries and shapes that 
are rather distinct from infinite, parallel and planar surfaces. Hence one 
faces the problem to compute the Casimir force between objects of general 
shape, arrangement and material decomposition. 

This article summarizes recent progress that has been proofed useful in 
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solving this problem for a variety of geometries. (For an overview of the 
development of related approaches, see Ref. 7.) In order to study Casimir 
forces in more general geometries, it turns out to be advantageous to de- 
scribe how fluctuating currents are induced on the objects by the scatter- 
ing of electromagnetic waves. This representation of the Casimir interaction 
was developed in Refs. 7-9. Each object is characterized by its on-shell elec- 
tromagnetic scattering amplitude. The separations and orientations of the 
objects are encoded in universal translation matrices, which describe how 
a solution to the source-free Maxwell's equations in the basis appropriate 
to one object looks when expanded in the basis appropriate to another. 
These matrices hence describe the electrodynamic interaction of the multi- 
pole moments associated with the currents and depend on the displacement 
and orientation of coordinate systems, but not on the shape and material 
of the objects themselves. The scattering amplitudes and translation ma- 
trices are then combined in a simple formula that allows efficient numerical 
and, in some cases, analytical calculations of Casimir forces and torques 
for a wide variety of geometries, materials, and external conditions. The 
approach applies to any finite number of arbitrarily shaped objects with 
arbitrary linear electromagnetic response at zero or finite temperature. 

To illustrate this general formulation, we provide some sample applica- 
tions, including results for the interaction between metallic objects for two 
spheres and for a sphere and a plane, taking into account the combined 
effect of shape and material properties at large distances. In addition, we 
provide examples for the non-additivity of the interaction by considering 
three objects (two spheres and a plane) and for the orientation dependence 
in the case of spheroids. The results are presented in form of analytical ex- 
pressions at large distances and as numerical results at smaller separations. 



2. Fluctuating currents and T-operators 

We consider the Casimir energy for neutral objects with electric and mag- 
netic susceptibilities. The partition function Z is defined through the path 
integral, which sums all configurations of the electromagnetic field (outside 
and inside the objects) with periodic boundary conditions in time between 
and T. The free energy F of the field at inverse temperature /3 is 

F{P) = - 1 p\ogZ{P). (1) 

The unrenormalizcd free energy generally depends on the ultraviolet cutoff, 
but cutoff-dependent contributions arise from the objects individually and 
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do not depend on their separations or orientations. Since we are only inter- 
ested in energy differences, we can remove these divergences by subtracting 
the energy of the system when the objects are in some reference configura- 
tion, see below. By replacing the time T by —ih/3, we obtain the partition 
function Z(f3) in 4D Euclidean space. In A° = gauge, the result is simply 
to replace the Matsubara frequencies oj n — by = icn n , where K n is 
the n th Matsubara frequency divided by c. The action is quadratic, so the 
modes with different n n decouple and the partition function decomposes 
into a product of partition functions for each mode. In the limit f3 — > oo, 
the sum J2 n >o turns mto an integral ^ f °° dn, and we have the ground 
state energy 



he 



/ dn log Z(k), (2) 
Jo 



with 



Z(k) = J VAVA* cxp -/3 J rfxE*(K,x) (u + ^V(k,x)^ E(k,x) 

(3) 

where we have used V x E = i^B to eliminate B in the action, and it 
is assumed that E is expressed by E = — c~ 1 d t A in terms of the vector 
potential A. This functional integral sums over configurations of A. This 
sum must be restricted by a choice of gauge, so that it does not include 
the infinitely redundant gauge orbits. We will choose to work in the gauge 
A = 0, although of course no physical results depend on this choice. Here 
we defined the Helmholtz operator 

i„(K)=I|^VxVx, (4) 

which is inverted by the Green's function that is defined by 

) =M< 3 )(x-x'). (5) 

The potential operator is 

V(k,x) = ln 2 (e(icK,x) - 1) + V x f — - 1 ) V x . (6) 

It is nonzero only at those points in space where the objects are located 
(e 7^ 1 or fi 1). At small frequencies, typical materials have e > 1 and 
/j, w 1, and V can be regarded as an attractive potential. 
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Next, we transform to a free field (with kernel Ho) by introducing fluctu- 
ating currents J that are confined to the objects. To perform this Hubbard- 
Stratonovich-like transformation we multiply and divide the partition func- 
tion of Eq. (3) by 



W= I 2?J2?J*| obj exp 



-p j dxJ*(x) •V" 1 (k,x)J(x) 



detV, (7) 



where | ob j indicates that the currents are defined only over the objects, i.e. 
the domain where V is nonzero (and therefore V -1 exists), and we have 
represented the local potential as a matrix in position space, V(k, x, x') = 
V(K,x)r5( 3 )(x — x'). We then change variables in the integration, J(x) = 
J'(x) + £V(k,x)E(x) and J*(x) = J'*(x) + iV(K,x)E'(x), to obtain 



Z(k) 



W 



VAVA* V3'V3'*\ , . x 

I obj 



exp 



-pjdxE*{ K ,x) (u {k) + ^V(/s,x)^ E(k,x) 



J'*(x) 



k,x)E*(k,x) V _1 (k,x) J'(x) 



(k, x)E(k, x) 



Wj 



VAVA* VJ'VJ' 



I obj 



exp 



dxE*H E + J'*V _1 J' + - (J'*E + J'E*) 



(8) 



Now the free electromagnetic field can be integrated out using H 1 = ac 2 Go, 
yielding 



Z(k) 



Zo 
W 



exp 



/ M'M'lobj 0) 
-P j dxdx' J'* (x) (G (k, x, x') + V" 1 (k, x)(5 3 (x - x')) J'(x') 



with Z a = JVAVA* exp[-/3 /dxE*H (K)E]. Both factors W and Z con- 
tain cutoff-dependent contributions but are independent of the separation 
of the objects. Hence these factors cancel and can be ignored when we con- 
sider a change in the energy due to a change of the object's separations with 
the shape and the material composition of the objects fixed. The kernel of 
the action in Eq. (9) is the inverse of the T-operator, i.e., T _1 = G + V -1 
which is equivalent to 



T = Vfll + GoV)- 1 . 



(10) 
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The Casimir energy at zero temperature (without the cutoff-dependent 
parts) is hence 

he r°° 

£ = — - / d/elogdetT. (11) 
27r Jo 

The determinant is here taken over the spatial indices x and x', which are 
restricted to the objects since T vanishes if x or x' are not on an object. 
To compute the determinant we start from the expression for T _1 which 
yields the reciprocal of the determinant. We decompose T -1 by introducing 
separate position space basis functions for each object. The projection of 
the currents onto this basis defines the object's multipole moments. This 
yields a division of T _1 into blocks where each block is labeled by an object. 

The off-diagonal blocks are given by Go only and describe the interaction 
of the multipoles on different objects. To see this we choose for each object 
individually an eigenfunction basis to expand the free Green's function, 

G («,x,x') =^E° ut («,x>)®E- g *(«;,x' < ) (12) 

a 

with regular solutions E^ cg and outgoing solutions E° ut of the free vector 
Helmholtz equation, where x< and x> denote the position with smaller 
and greater value of the "radial" variable of the separable coordinates. The 
multipole moments of object j are then Qj ta (n) = J rfxJj ;(k, x)E^ cg * (k, x). 
Regular solutions form a complete set and hence outgoing solutions can 
be expanded in terms of regular solutions except in a region (enclosed by a 
surface of constant radial variable) that contains the origin of the coordinate 
system of object i. This expansion defines the translation matrices Ug Q via 

E° ut (« )Xi ) = (k.X^E^k.x,-) , (13) 



where the definition of the coordinates is shown in Fig. 1. The free Green's 
function then becomes 

Go^x^'J^E^^x^U^KX^E^'^x;) (14) 

a, 13 

so that the off-diagonal blocks of T _1 are given by the translation matrices. 
Equivalent translation matrices can be defined between two sets of regular 
solutions as is necessary for one object inside another, see Ref. 7. 

The diagonal blocks of T" 1 are given by the matrix elements of the 
T-operators Tj of the individual objects. By multiplying T -1 by the T- 
operator without the off-diagonal blocks which can interpreted as de- 
scribing a reference configuration with infinite separations between the ob- 
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Fig. 1. Geometry of the configuration. The dotted lines show surfaces separating the 
objects on which the radial variable is constant. The translation vector = x^ — Xj = 
— Xji describes the relative positions of the two origins. 

jects, one finds that (for objects outside each other) the diagonal blocks 
are given by the inverse of the matrix representing Tj in the basis E" s . 7 
The physical meaning of this matrix follows from the Lippmann-Schwinger 
equation for the full scattering solution E Q (re, x), 

E a («,x) = E^( K ,x)-G V,E Q ( K ,x) = E- g ( K ,x)-G T J E-s( K ,x) . (15) 

Using the expansion of Eq. (12), the solution sufficiently far away from the 
object (i.e., for positions that have a radial variable larger than any point 
on the object) can be expressed as 

E a («,x) = E r a es («,x) - ]TE^( K ,x) f E^(«,xOT,( K )E- g ( K ,x')dx', 

(16) 

where the integral defines the scattering amplitude Wj t p a (K,) of object j. It 
can be obtained, e.g., from matching boundary conditions at the surface of 
a dielectric object. 

The Casimir energy (without cutoff-dependent contributions from W 
and Z ) can now be expressed as 

he f°° 

£= 2^Jo dKl °g det ( MM -)> ( 17 ) 

where 

'F^ 1 U 12 U 13 • 

U 21 F^ 1 U 23 • • • | (18) 
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and M^ 1 is the block diagonal matrix diag(Fi, F2, ■ ■ ■ )• For ^ ne case °f two 
objects this expressions simplifies to 

£ = — / d/dogdet (I - FiU 12 F 2 U 21 ) . (19) 

In order to obtain the free energy at nonzero temperature instead of the 
ground state energy, we do not take the limit (5 — > 00 in Eq. (I). 3 Instead, 
the integral ^ J °° dn is replaced everywhere by ^J2 n , where cn n = 
with n = 0,1,2,3... is the nth Matsubara frequency. A careful analysis 
of the derivation shows that the zero frequency mode is weighted by 1/2 
compared to the rest of the terms in the sum; this modification of the sum 
is denoted by a prime on the summation symbol. 

3. Applications 

In this section we demonstrate the applicability of the method through some 
examples. Due to the lack of space, we only present the final analytical and 
numerical results that all follow from Eq. (17) or Eq. (19) by truncation of 
the matrices at some order of partial waves, i.e., by considering only a finite 
set of basis functions. At asymptotically large distances, the interaction 
only depends on the dipole contribution while with drecreasing distance 
the number of partial waves has to be increased. Below we will provide 
results both in form of a asymptotic series in the inverse separation and 
numerical results for a wide range of distances. 

3.1. Sphere-plane 

First, we consider the sphere-plate geometry that has been employed in 
the majority of recent experiments. At large distances, the energy can be 
expanded in an asymptotic series in the inverse separation. For a dielectric 
sphere in front of perfectly reflecting mirror with sphere-center to mirror 
separation L the Casimir energy is 

* = — a ^j4 + - a * + 2 7i3 - 27i M 3 )^ 

+ T(r24 [23K)2 _ 14<af + 23K)2 + 2160(7 " 4 _ 7 " )] h (20) 

+ 72^0 [572K " a * ] + 675 (9(7 " 5 " 7l " } " 55(723 " 723))] lb + ■ • ' } ' 

where af, af are the static electric and magentic multipole polarizabilities 
of the sphere of order I (I = 2 for dipoles), and the coefficients jf n , ~ff n 
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describe finite- frequency corrections to these polarizabilities, i.e., terms ~ 
K 2i+n m ^ e Iow-k expansion of the T-matrix clement for the I th partial 
wave. Notice that the first three terms of the contribution at order L~ 7 
have precisely the structure of the Casimir- Polder interaction between two 
atoms with static dipole polarizabilities af and ctf but it is reduced by a 
factor of 1/2 8 . This factor and the distance dependence ~ L of this term 
suggests that it arises from the interaction of the dipole fluctuations inside 
the sphere with those inside its image at a distance 2L. The additional 
coefficient of 1/2 in the reduction factor (1/2) (1/2 7 ) can be traced back to 
the fact that the forces involved in bringing the dipole in from infinity act 
only on the dipole and not on its image. If the sphere is also assumed to be 
a perfect reflector, the energy becomes 

£ = -tYM t - ( 21 ) 



where the coefficients up to order 1/L 11 are 

i 9 25 3023 

Od = , Or; = (J, Or = , 07 = 

4 16' 5 ' 6 32' 7 4096 
12551 , 1282293 



9600 ' ° 163840 ' 
32027856257 _ 39492614653 

10 ~ 722534400 ' 11 ~ 412876800 ' ^ ' 

Our method can be also employed to study the material dependence of 
the interaction. When the sphere and the mirror are described by a simple 
plasma model, we can obtain the interaction energy again from Eq. (19) by 
substituting the dielectric function on the imaginary frequency axis, 

= 8 ■ (23) 

into the T-matrices of sphere and mirror. From this we get at large sepa- 
rations 



he 
ir 



d3 n>4 

U{\ P /R)^ + MXp/R)^ + o(L- e ) 



with the functions 



, . , 9 9 o 9 , 2tt 

/4(Z)= 16 + 64^ Z -32^ C ° th T 
, . . 13 21 , 21 2 , 2tt 

/5W = -20^-80t^ + 40^ Z C ° th T 



(24) 



(25) 
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It is interesting that the amplitude fa of the leading term is not universal 
but depends on the plasma wavelength X p . Only in the two limits X p /R — > 
and Xp/R — > oo the amplitude assumes material independent values, 9/16 
and 3/8, respectively. The first limit describes perfect reflection of electric 
and magnetic fields at arbitrarily low frequencies and hence agrees with the 
result of Eq. (21). The change to the second amplitude for large X p can be 
understood when one considers a London superconductor that is described 
at zero temperature by the plasma dielectric function. 10 If one associates 
X p with the penetration depth, the perfect reflector limit results from the 
absence of any field penetration while the second limit corresponds to a 
large penetration depth and hence the suppression of the magnetic mode 
contribution to the Casimir energy, explaining the reduced amplitude of 
3/8. The latter result follows also when the objects are considered to be 
normal metals, described by the Drude model dielectric function 

e v {icK) = l+ ff* . (26) 

From this function we get for a sphere and a mirror made of a Drude metal 
the asymptotic energy 

_ he h R 3 77 R 3 ( c O(L-t-) 

(27) 

In fact, one observes that the leading term is universal and agrees with the 
Xp — > oo limit of the plasma model. Note that the result of Eq. (27) does not 
apply to arbitrarily large dc conductivity a. The conditions for the validity 
of Eq. (27) can be written as L > R, L > c/a and L > aR 2 /c. The 
above results demonstrate strong correlations between shape and material 
since for two parallel, infinite plates, both the plasma and the Drude model 
yield at large separations the same (universal) result as a perfect mirror 
description. 

In order to study short separations, Eq. (19) has to be evaluated nu- 
merically by including sufficiently many partial waves. The result of an 
extrapolation from I — 29 partial waves is shown in Fig. 2 in the perfect 
reflection limit. 11 At small separations the result can be fitted to a power 
law of the form 

2 



£ — £pFA 



(28) 



with £pfa and d defined in Fig. 2. The coefficients 0j measure corrections 
to the proximity force approximation and are obtained from a fit of the 
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function of Eq. (28) to the data points for the four smallest studied sepa- 
rations. We find 6»i = -1.42 ± 0.02 and 9 2 = 2.39 ± 0.14. This result is in 
agreement with numerical findings in Ref. 12 but is in disagreement with 
an asymptotic expansion for small distances. 13 The latter yields 6\ = —5.2 
and very small logarithmic corrections that however can be ignored at the 
distances considered here. The origin of this discrepancy is currently un- 
clear but might be related to the applicability of the asymptotic expansion 
to only much smaller distances than accessible by current numerics. 




February 28, 2010 16:48 WSPC - Proceedings Trim Size: 9in x 6in qfext09 



11 




o - 



E |H/L->0 M 



cancel 

1~- 



add 



I H/L^ 



1 |i 



0* 



'.(90 



Fig. 3. Left: Geometry of the two-sphere/atom and sidewall system. Shown are also 
the mirror images (grey) and two- and three-body contributions (solid and dashed curly 
lines, respectively). Right: Typical orientations of electric (E) and magnetic (M) dipoles 
and image dipoles for H/L — > and H/L — > co. 



3.2. Three-body effects 

Casimir interactions are not pair-wise additive. To study the consequences 
of this property, we consider the case of two identical, general polarizable 
objects near a perfectly reflecting wall in the dipole approximation, see 
Fig. 3. This situation applies to ground state atoms and also to general 
objects at large separations. The separation between the objects is L and 
the separation of each of them from the wall is H. In dipole approximation, 
the retarded limit of the interaction is described by the static electric (a z , 
ay) and magnetic (/3 Z , /3m) dipole polarizabilities of the objects which can 
be different in the directions perpendicular (z) and parallel (||) to the wall. 
In the absence of the wall the potential for the two polarizable objects is 
given by the well-known Casimir-Polder (CP) potential 

he r 1 
^2,|(i) = -^^[33ajj +13a z -14q||/3 z + (ao/3)J , (29) 

The L-dependent part of the interaction energy in the presence of the wall 
is 

S 9 o(L,H)=S 2A (L)+S 2A (D,L) + £ 3 (D,L) (30) 

with D = V L 2 + AH 2 . The change in the relative orientation of the objects 
with I = L/D leads to the modified 2-body CP potential 

he r 

£ 2 \{D,L) = -— — = 26a 2 +20a 2 z -U£ 2 (4a 2 - 9ana z + 5a 2 z ) (31) 

OTTjJ'l 11 11 

+ 63f(a|| - a z) 2 - 14(a||^||(I-^ 2 ) +£ 2 a n /3 z ) + (ao/?)] . 
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The 3-body energy £^(D,L) describes the collective interaction between 
the two objects and one image object. It is given by 

£ 3 (D, L) = jfe ^— L— ^ [(3^ + 15£ 5 + 28£ 4 + 20£ 3 + U 2 -U-l) 
x [a\ - /3f ) - (3£ 6 + 15^ + 24£ 4 - 10£ 2 -5t-l) [a\ - f3 2 z ) 
+ 4(f 4 + 5£ 3 +f 2 ) (a t 0\\-a\\P x )]. (32) 

It is instructive to consider the two limits H -C L and H ^> L. For H <C L 
£oo_ turns out to be the CP potential of Eq. (29) with the replacements a z — > 
2a z , au — » 0, f3 z —¥ 0, -4- 2/3m. The 2-body and 3-body contributions add 
constructively or destructively, depending on the relative orientation of a 
dipole and its image which together form a dipolc of zero or twice the 
original strength (see Fig. 3). 

For H » L the leading correction to the CP potential of Eq. (29) comes 
from the 3-body energy The energy then becomes (up to order H~ 6 ) 



he 

£ 00 (L,H)=£ 2il (L) + - 

7T 



+ " or^ (« P) 



4:L 3 H 4 8LH e 



■ (33) 



The signs of the polarizabilitics in the leading term <~ H~ 4 can be under- 
stood from the relative orientation of the dipole of one atom and the image 
dipole of the other atom, see Fig. 3. If these two electric (magnetic) dipoles 
are almost perpendicular to their distance vector they contribute attrac- 
tively (repulsively) to the potential between the two original objects. If 
these electric (magnetic) dipoles are almost parallel to their distance vector 
they yield a repulsive (attractive) contribution. For isotropic polarizabili- 
ties the leading term of Eq. (33) vanishes and the electric (magnetic) part 
~ H~ 6 of the 3-body energy is always repulsive (attractive). 

Next, we study the same geometry as before but with the objects as- 
sumed to be two perfectly reflecting spheres of radius R. The lengths L and 
H are measured now from the centers of the spheres, see Fig. 3. Here we do 
not limit the analysis to large separations but consider arbitrary distances 
and include higher order multipole moments than just dipole polarizability. 
For R<^ L, H and arbitrary H/L the result for the force can be written as 



3=6 



The functions fj can be computed exactly. We have obtained them up to 
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(35) 



(36) 



For H > L one has f 6 (h) = -1001/16 + 3/(4h 6 ) + C(^ 8 ), h{h) = 
-71523/160 + 39/(80/i 6 ) + <D{h~ 8 ) so that the wall induces weak repul- 
sive corrections. For H < L, f 6 (h) = -791/8 + 6741h 2 /8 + 0(h 4 ), 
fs(h) = -60939/80 + 582879/i 2 /80 + 0{h 4 ) so that the force amplitude 
decreases when the spheres are moved a small distance away from the wall. 
This proves the existence of a minimum in the force amplitude as a function 
of H/R for fixed, sufficiently small R/L. We note that all fj(h) are finite 
for h — > oo but some diverge for Ii4 0, e.g., fg <~ /n <~ h~ 3 , making them 
important for small H. 

To obtain the interaction at smaller separations or larger radius, we have 
computed the energy £oo and force F = —d£oo_/dL between the spheres 
numerically. 14 In order to show the effect of the wall, we plot the energy and 
force normalized to the results for two spheres without a wall. Fig. 4 shows 
the force between the two spheres as a function of the wall distance for fixed 
L. When the spheres approach the wall, the force first decreases slightly if 
R/L < 0.3 and then increases strongly under a further reduction of H. 
For R/L > 0.3 the force increases monotonically as the spheres approach 
the wall. This agrees with the prediction of the large distance expansion. 
The expansion of Eq. (34) with j = 10 terms is also shown in Fig. 4 for 
R/L < 0.2. Its validity is limited to large L/R and not too small H/R; it 
fails completely for R/L > 0.2 and hence is not shown in this range. 

3.3. Orientation dependence 

In this section we investigate the shape and orientation dependence of the 
Casimir force using Eq. (19). As examples we focus on ellipsoids, comput- 
ing the orientation dependent force between two spheroids, and between a 
spheroid and a plane. 15 For two anisotropic objects, the CP potential of 
Eq. (29) must be generalized. In terms of the Cartesian components of the 
standard electric (magnetic) polarizability matrix a (/3), the asymptotic 



February 28, 2010 16:48 WSPC - Proceedings Trim Size: 9in x 6in qfext09 



14 




1 2 (if - R)/R 3 4 5 



Fig. 4. Numerical results for the force (dots) between two spheres as function of the 
sidewall separation H/R for different sphere separations R/L. Shown are also the ana- 
lytical results of Eq. (34), including terms up to j = 10 for R/L < 0.2 (solid curves). 
Inset: Magnification of the nonmonotonicity. 



large distance potential of two objects (with the z axis pointing from one 
object to the other), can be written as 

F — — — — il ( 12,12 , o i 2 \ 
^ ^7 g^. | y^-xx^xx ' ^yy^yy ' ^ a xy^xy) 

+ 20 a\ z a 2 zz - 30 {a\ z al z + a\ z a\ z ) + (a -> /3) (37) 

— ^ {pt-xxfiyy + a yyPxx ~ ^ a xyPxy) + (1 ^ 2) 



For the case of an ellipsoidal object with static electric permittivity e and 
magnetic permeability /j,, the polarizability tensors are diagonal in a basis 
oriented to its principal axes, with elements (for i € {1,2,3}) 

n .o = y e-l 0= V m-1 m] 

u 4tt 1 + (e - l)n* ' Pvl 4tt 1 + (/i - l)n; ' k ; 

where F = 47rr 1 r2r 3 /3 is the ellipsoid's volume. In the case of spheroids, 
for which r\ = T2 = R and r% = L/2, the so-called depolarizing factors can 
be expressed in terms of elementary functions, 

1 _ n 3 1 - e 2 /. 1 + e \ 

rn=n 2 = ^— , ^ = ^3- I log— -2e], (39) 
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where the eccentricity e = yj 1 — is real for a prolate spheroid (L > 2R) 
and imaginary for an oblate spheroid (L < 2R). The polarizability tensors 
for an arbitrary orientation are then obtained as a = lZ~ l a°1Z, where 1Z is 
the matrix that rotates the principal axis of the spheroid to the Cartesian 
basis, i.e. 72.(1,2,3) — > (x,y,z). Note that for rarefied media with e ~ 1, 
fi ~ 1 the polarizabilities are isotropic and proportional to the volume. 
Hence, to leading order in e — 1 the interaction is orientation independent 
at asymptotically large separations, as we would expect, since pairwise sum- 
mation is valid for e — 1 <C 1. In the following we focus on the interesting op- 
posite limit of two identical perfectly reflecting spheroids. We first consider 
prolate spheroids with L » R. The orientation of each "needle" relative to 
the line joining them (the initial z-axis) is parameterized by the two angles 
(9,ip), as depicted in Fig. 5. Then the energy is 



he 



5L 6 



d7 I 1152tt (In I -if 



cos 2 9\ cos 2 9 2 



(40) 



13 3 
+ — cos 2 ip sin 2 6i sin 2 9 2 — - cos tp sin 29 1 sin 29 2 
20 8 



O 



L 4 R 2 



where ip = ipi — tp 2 - It is minimized for two needles aligned parallel to 
their separation vector. At almost all orientations the energy scales as L 6 , 
and vanishes logarithmically slowly as R — > 0. The latter scaling changes 
when one needle is orthogonal to z (i.e. 9\ = tt/2), while the other is either 
parallel to z (9 2 — 0) or has an arbitrary 9 2 but differs by an angle 7r/2 in 
its rotation about the z-axis (i.e. ^p\—ip2 = 7r /2). In these cases the energy 
comes from the next order term in Eq. (40), and takes the form 



he 



L^R 2 



1152nd 7 In £ 

IX 



(73 + 7 cos 26*2 ) 



(41) 



which shows that the least favorable configuration corresponds to two nee- 
dles orthogonal to each other and to the line joining them. 

For perfectly reflecting oblate spheroids with R ^> L/2, the orientation 
of each "pancake" is again described by a pair of angles (9, tp), as depicted 
in Fig. 6. To leading order at large separations, the energy is given by 



hc( R 6 
^\144tt 3 



765 - 5(cos 20i + cos 29 2 ) + 237 cos 29 x cos 26> 2 



+ 372 cos 2ip sin 2 9 1 sin 2 6 2 - 300 cos tp sin 29 1 sin 26> 2 



+ 0(R 5 L) 



(42) 
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The leading dependence is proportional to i? 6 , and does not disappear for 
any choice of orientations. Furthermore, this dependence remains even as 
the thickness of the pancake is taken to zero (L — > 0). This is very different 
from the case of the needles, where the interaction energy vanishes with 
thickness as In - 1 (L/R). The lack of L dependence is due to the assumed 
perfectly reflectivity. The energy is minimal for two pancakes lying on the 
same plane {Q x = 9 2 = tt/2, ip = 0) and has energy -he (173/18tt 3 )R 6 /d 7 . 
When the two pancakes are stacked on top of each other, the energy is 
increased to — he (62/9ir 3 )R e /d 7 . The least favorable configuration is when 
the pancakes lie in perpendicular planes, i.e., 9i = n/2, 6 2 = 0, with an 
energy -he (ll/3ir 3 )R 6 /d 7 . 




Fig. 5. (Color online) Orientation of a prolate (cigar-shaped) spheroid: The symmetry 
axis (initially the z-axis) is rotated by 9 about the rr-axis and then by ip about the z-axis. 
For two such spheroids, the energy at large distances is give by Eq. (40). The latter is 
depicted at fixed distance d. and for ipi = -02, by a contour plot as function of the angles 
01, 82 for the x-axis rotations . Minima (maxima) are marked by filled (open) dots. 



For an anisotropic object interacting with a perfectly reflecting mirror, 
at leading order the CP potential generalizes to 

he 1 

£ = -^-tr(a-(3)+<D(d- 5 ), (43) 

which is clearly independent of orientation. Orientation dependence in this 
system thus comes from higher multipoles. The next order also vanishes, 
so the leading term is the contribution from the partial waves with I = 3 
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Fig. 6. (Color online) As in Fig. 5 for oblate (pancake-shaped) spheroids, with a contour 
plot of energy at large separations. 



for which the scattering matrix is not known analytically. However, we 
can obtain the preferred orientation by considering a distorted sphere in 
which the radius R is deformed to R + Sf("&,<p). The function / can be 
expanded into spherical harmonics Ki m ($, tp), and spheroidal symmetry can 
be mimicked by choosing / = lao($) (f). The leading orientation dependent 
part of the energy is then obtained as 

^-""S^^ ""<*>■ (44) 

A prolate spheroid (6 > 0) thus minimizes its energy by pointing towards 
the mirror, while an oblate spheroid (5 < 0) prefers to lie in a plane per- 
pendicular to the mirror. (We assume that the perturbative results are not 
changed for large distortions.) These configurations are also preferred at 
small distances d, since (at fixed distance to the center) the object reori- 
ents to minimize the closest separation. Interestingly, the latter conclusion 
is not generally true. In Ref. 15 it has been shown that there can be a 
transition in preferred orientation as a function of d in the simpler case of 
a scalar field with Neumann boundary conditions. The separation at which 
this transition occurs varies with the spheroid's eccentricity. 



3.4. Material dependence 

In this section we shall discuss some characteristic effects of the Casimir 
interaction between metallic nano-particles by studying two spheres with 
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finite conductivity in the limit where their radius R is much smaller than 
their separation d. We assume further that R is large compared to the 
inverse Fermi wave vector ir/kp of the metal. Since typically ir/kp is of 
the order of a few Angstrom, this assumption is reasonable even for nano- 
particlcs. Theories for the optical properties of small metallic particles 16 
suggest a Drude dielectric function 

-, , v{icn) 
e(icn) = 1 + 4tt — , (45) 

CK 

where a (ion) is the conductivity which approaches for k — > the dc con- 
ductivity Ode- For bulk metals Ode — w^t/Ait where lo p = Ae 2 k F /'iiim e 
is the plasma frequency with electron charge e and electron mass m e , and 
r is the relaxation time. With decreasing dimension of the particle, Udc{R) 
is reduced compared to its bulk value due to finite size effects and hence 
becomes a function of R. 16 In analogy to the result for a sphere and a plate 
that are described by the Drude model, we obtain for the large distance 
expansion of the energy the result 

23 R e f Ra dc (R) 45 c \ R 7 

As in the sphere-plate case, the leading term is material independent but 
different from that of the perfect metal limit (where the amplitude is 
143/1671") since only the electric polarization contributes. At next order, 
the first and second terms in the parentheses come from magnetic and elec- 
tric dipolc fluctuations, respectively. The term ~ 1/L S is absent in the 
interaction between perfectly conducting spheres. The limit of perfect con- 
ductivity, Ode — > oo cannot be taken in Eq. (46) since this limit does not 
commute with the large L expansion. 

In order to estimate the effect of finite conductivity and its dependence 
on the size of the nano-particle, we have to employ a theory that can de- 
scribe the evolution of <7dc(R) with the particle size. A theory for the di- 
electric function of a cubical metallic particle of dimensions R ^> n/kp 
has been developed within the random phase approximation in the limit of 
low frequencies <C c/R. 16 In this theory it is further assumed that the dis- 
creteness of the electronic energy levels, and not the inhomogeneity of the 
charge distribution, is important. This implies that the particle responds 
only at the wave vector of the incident field which is a rather common 
approximation for small particles. From an electron number-conserving re- 
laxation time approximation the complex dielectric function is obtained 
which yields the size-dependent dc conductivity for a cubic particle of vol- 
ume o 3 . 16 It has been shown that the detailed shape of the particle does 
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Fig. 7. Dimensionlcss dc conductivity a^ c {R) in units of e 2 /2har> (with Bohr radius 

ao) for a Aluminum sphere with ep = 11.63eV, ir/kp = 1.8A and r = 0.8 • 10~ 14 sec as 
function of the radius R, measured in units of ir/kp. Also shown is the corresponding ratio 
RiJdc (R) l c that determines the Casimir interaction of Eq. (46) . The bulk dc conductivity 
(Td c (oo) = 17.66 is indicated by the dashed line. 

not matter much, and we can set a — (47r/3) 1 / 3 i? which defines the volume 
equivalent sphere radius R. For it /hp — a the nano particle ceases to be 
conducting, corresponding to a metal-insulator transition due to the local- 
isation of electrons for particles with a size of the order of the mean free 
path. It is instructive to consider the size dependence of Od c {R) and of the 
Casimir interaction for a particular choice of material. Following Rcf. 16, 
we focus on small Aluminum spheres with Fermi energy ep = 11.63eV 

o 

and t = 0.8 • 10 sec. These parameters correspond to Tr/k F = 1.8 A 
and a plasma wavelength X p = 79nm. It is useful to introduce the di- 
mensionless conductivity &d c (R), which is measured in units of e 2 /2ha 
with Bohr radius do, so that the important quantity of Eq. (46) can be 
written as R<Jd c (R)/c — (a/2)(R/a )<Jd c (R) where a is the fine-structure 
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constant. The result is shown in Fig. 7. For example, for a sphere of radius 
B = lOnm, the dc conductivity is reduced by a factor w 0.15 compared to 
the bulk Drude value. If the radius of the sphere is equal to the plasma 
wavelength X p , the reduction factor w 0.8. These results show that shape 
and material properties are important for the Casimir interaction between 
nano-particles. Potential applications include the interaction between dilute 
suspensions of metallic nano-particles. 

3.5. Further extensions 

The general result of Eq. (17) and its extensions described in Ref. 7 have 
been recently applied to a number of new geometries and further appli- 
cations are under way. Examples include so-called interior configurations 
with an object contained within an otherwise empty, perfectly conducting 
spherical shell. 17 For this geometry the forces and torques on a dielectric or 
conducting object, well separated from the cavity walls, have been deter- 
mined. Corrections to the proximity force approximation for this interior 
problem have been obtained by computing the interaction energy of a finite- 
size metal sphere with the cavity walls when the separation between their 
surfaces tends to zero. Eq. (17), evaluated in parabolic cylinder coordinates, 
has been used to obtain the interaction energy of a parabolic cylinder and 
an infinite plate (both perfect mirrors), as a function of their separation 
and inclination, and the cylinder's parabolic radius. 18 By taking the limit 
of vanishing radius, corresponding to a semi-infinite plate, the effect of edge 
and inclination could be studied. 
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